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ABSTRACT 


This  goal  of  this  research  is  to  characterize  two  micro  electro  mechanical 
system  (MEMS)  based  directional  sound  sensors  with  solid  and  perforated 
wings.  The  design  of  the  sensors  was  based  on  the  structure  of  Ormia  ochracea 
fly's  hearing  system  that  has  highly  directional  hearing  through  mechanical 
coupling  of  the  eardrums.  The  sensors  are  made  of  10  micron  thick  single  crystal 
silicon  layer  with  dimensions  1  x  2  mm2.  The  sensors  were  fabricated  using 
SOIMUMPs  process  available  through  MEMSCAP  foundry  service.  The 
characteristics  of  the  two  sensors  were  simulated  COMSOL  finite  element 
software  and  responses  to  incident  sound  at  different  angles  were  measured 
using  a  laser  vibrometer.  Both  sensors  showed  good  sound  coupling  and 
measured  and  simulated  frequency  responses  are  in  good  agreement.  The 
sensor  with  perforated  wings  was  found  to  have  faster  response  compared  with 
that  of  the  solid  wings  primarily  due  to  lower  mass  and  higher  damping.  The 
measurements  showed  good  sensitivity  to  the  direction  of  sound  as  predicted 
from  the  modeling. 
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I.  INTRODUCTION 


A.  INTRODUCTION 

The  new  design  of  generation  4th  of  directional  microphone  is  designed  to 
make  an  improvement  from  the  previous  efforts  of  understanding  the  sensitivity 
of  directional  hearing.  Many  small  animals  depend  on  the  ability  to  localize  the 
sound  sources  (Richard  R.  Fay  and  Arthur  N.  Popper,  1992).  All  animals  having 
two  tympanal  ear  (eardrums)  to  localize  sound  by  processing  interaural 
differences  either  in  the  time  of  arrival  or  level  of  the  acoustic  pressure.  When 
the  size  of  the  animal  is  very  small  to  relative  to  the  sound  wavelength,  these 
interaural  differences  can  be  too  small  to  permit  accurate  processing  by  the 
central  nervous  system  to  allow  sound  source  localization.  By  suitable  coupling 
of  the  motions  of  the  tympana,  it  is  possible  for  each  to  respond  preferentially  to 
sound  from  certain  directions.  This  paper  is  to  analysis  of  micromechanical 
device  inspired  by  the  coupled  ears  of  the  parasitoid  fly,  Ormia  ochracea,  and  its 
usefulness  as  a  has  inspired  a  novel,  approach  to  constructing  small  device  and, 
directional  sound  receivers  (Daniel  Robert  and  Ronald  R.  Hoy,  1998).  The 
objective  of  this  directional  microphone  is  to  design  to  make  an  improvement 
from  previous  efforts  (Antonios  Dritsas,  June  2008)  to  understand  the  sensitivity 
of  this  method  of  directional  hearing. 

This  paper  is  a  continuation  development  from  the  previous  model,  which 
is  trying  to  improve  the  effect  for  acoustic  directional  sensing.  The  generation  4th 
of  Silicon  devices  studied  shows  significant  potential  for  enabling  the 
development  of  novel  sensors  for  sound  and  vibration  and  has  focused  on  the 
fabrication  of  small  design  for  sound.  This  experiment  is  to  demonstrate  the 
improvement  in  speech  intelligibility  in  noise  to  achieve  directional  acoustic 
response  and  I  have  seen  very  good  improvements  in  solid  state  in  its 
performance  in  terms  of  quality  factor  and  the  results  presented  here  will  also 
help  enable  the  future  development  of  a  good  design  for  direction  finding. 
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B.  BACKGROUND 


Some  large  insects  creatures  can  use  only  sound  amplitude  pressure 
difference  and  arrival  time  difference  for  their  directional  hearing  and  others  such 
as  insects  can  make  use  of  both  amplitude  and  phase.  In  humans  head 
(diameter  «  17  cm),  the  difference  in  the  time  of  arrival  at  the  two  ears  is 
approximately  0.5  ms.  In  insects  the  part  of  the  body  carrying  the  ears  is  about 
10  to  50  times  smaller  than  the  human  head,  and  large  differences  in  sound 
pressure  at  the  ears  only  exist  at  high  frequencies.  Most  insects  thus  have  great 
difficulty  in  determining  the  direction  of  sound  incidence  by  measuring  the 
difference  in  time  when  the  sound  waves  reach  both  the  left  and  right  ears. 
Alternative  ways  without  involving  the  pressure  component  of  sound  is  by 
studying  the  Mechanical  analogue  of  the  tympanal  system. 

C.  THE  PROBLEM  OF  SMALL  SIZE  FOR  DIRECTIONAL  HEARING 

Animals  use  two  basic  acoustics  cues,  interaural  time  and  amplitude 
differences,  for  the  directional  binaural  detection  of  a  sound  source.  From  the 
physical  laws  of  sound  propagation  in  air  and  the  diffraction  of  sound  around 
solid  bodies,  it  is  apparent  that  very  small  animals  such  as  insects  face 
limitations  in  their  ability  to  extract  detectable  interaural  time  and  amplitude 
acoustic  cues  from  the  incident  sound  field  (Middlebrooks  and  Green,  1991).  In 
theory,  given  the  small  body  size  of  the  ormiine  flies  and  the  very  short  interaural 
distance,  the  generation  of  reliable  interaural  acoustic  cues  seems  highly 
problematic. 

It  is  important  to  mention,  in  passing,  that  time  and  amplitude  acoustic 

cues  are  not  necessarily  required  to  localize  an  incident  sound.  Monaural  sound 

localization  by  means  of  spectral  cues  is  possible  in  humans  under  certain 

circumstances.  The  processing  of  spectral  information  for  directional  hearing  in 

insects  has  not  been  given  much  attention  and  can  therefore  not  be  excluded 

from  consideration.  As  the  speed  of  sound  344  m/s  in  air,  an  interaural  distance 

of  1  cm  will  only  generate  an  interaural  time  difference  of  the  incident  sound 
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wave  (ITD)  of  about  30  ps.  Such,  and  even  shorter,  time  delays  admittedly  pose 
a  severe  challenge  to  the  nervous  system  in  its  typical  temporal  range  of 
operation  on  the  order  of  a  millisecond.  The  other  basic  cue,  the  interaural 
intensity  differences,  has  operational  size  limits  that  are  related  to  the  wavelength 
of  the  incident  sound  and  the  size,  shape,  and  density  of  the  body  or  head 
carrying  the  ears.  It  is  generally  accepted  from  acoustical  diffraction  theory  that  a 
ratio  of  1:10  between  size  and  wavelength  does  not  give  rise  to  measurable 
diffractive  effects  (Morse  and  Ingard,  1968). 

For  example,  a  spherical  body  of  7  mm  in  diameter  does  not  significantly 
diffract  sounds  below  a  frequency  of  5  kHz  (with  a  wavelength  of  68  mm),  and 
therefore  very  little  or  even  no  difference  in  interaural  sound  pressure  occurs.  It 
uses  5  kHz  as  an  example  because  it  is  the  dominant  frequency  of  field  crickets’ 
calling  song  and  therefore  is  salient  to  its  acoustic  parasitoid  as  well. 

In  theory,  given  the  small  body  size  of  the  ormiine  flies  (a  few  millimeters 
in  breadth)  and  the  very  short  interaural  distance  (1.68  mm  total  width  for  both 
tympanal  membranes  taken  together),  the  generation  of  reliable  (and 
physiologically  relevant)  interaural  acoustic  cues  seems  highly  problematic. 
Practically,  how  does  the  fly  acoustically  localize  her  cricket  host  singing  at  5 
kHz?  In  the  fly’s  case,  the  body  size  to  wavelength  ratio  is  about  1:130,  thus 
precluding  the  effects  of  diffraction  as  the  source  of  interaural  intensity  difference. 
An  angle  of  incidence  of  90  degree  relative  to  the  normal  is  the  best  case  for 
sound  localization,  and  sound  will  travel  from  one  side  to  the  other  of  the  most 
lateral  margins  of  the  tympanal  membranes  (1.68  mm  apart)  in  about  4.9  ps.  In 
sensory  organs  is  about  1.5  ps,  given  their  separations  of  520  ps.  Actual 
measurements  made  with  custom-made  probe  microphones  confirm  these 
theoretical  arguments.  At  5  kHz,  the  interaural  intensity  difference  is  too  small  to 
be  measured,  (<1dB)  either  across  the  fly’s  body  or  at  the  tympanal  membranes. 
The  maximal  interaural  time  delay  measured  by  two  phase-calibrated  probe 
microphones  positioned  directly  in  front  of  the  tympanal  membranes  is  1 .45  ps. 
[Robert,  Amoroso,  and  Hoy,  1992], 
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II.  THE  DESCRIPTION  OF  ORMIINE  FLIES 


The  ears  of  ormiine  flies  are  positioned  on  the  ventral  prothroax,  above 
and  anterior  to  the  first  pair  of  legs  and  just  behind  the  head  (Figure  1  .A)  (Lakes- 
Harlan  and  Heller  1992;  Robert,  Amoroso,  and  Hoy  1992).  The  general  design 
feature  of  the  ormiine  hearing  organs  is  consistent  with  the  design  of  other 
insectan  tympanal  organs  (Robert,  Read,  and  Hoy  1994;  Hoy  and  Robert  1996). 
A  tympanal  ear  in  insects  essentially  consists  of  three  basic,  major  morphological 
components: 

1.  A  localized  thinning  of  the  ‘cuticle’  to  provide  a  tympanal 
membrane,  which  is  formed  by  the  apposition  of  a  modified,  thin 
exocuticle  and  a  tracheal  air  sac. 

2.  An  air  chamber  backing  the  tympanal  membrane,  which  is  formed 
by  the  associated  tracheal  air  sac. 

3.  A  mechanoreceptive  sensory  organ  of  the  scolopidial  type,  which  is 
indirect  contact  with  either  the  tympanal  membrane  or  the  modified 
tracheal  system  abutting  it. 


Figure  1 .  External  auditory  anatomy  of  Ormia  ochracea.  In  ormiine  flies,  the 
tympanal  ears  are  located  between  the  first  pair  of  legs  and  the  base  of 
the  neck.  The  light  scanning  micrographs  show  a  semilateral  view  of  an 
intact  fly  (A),  and  (B)  of  a  fly  with  the  head  removed  to  see  the  prostenal 
hearing  organs.  Co,  prothoracic  coax;  N,  neck;  PTM,  prosternal  tympanal 
membranes;  Pb,  probasisternum;  MSP,  mesothoracic  spiracle  (Lakes 
Harlan  and  Heller  1992;  Robert,  Amoroso,  and  Hoy  1992). 
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A  series  of  eight  morphological  specializations  that  are  particular  to 
ormiines  have  been  identified: 

1.  Enlargement  of  the  prosternal  membranes  (PM),  giving  rise  to  thin 
prosternal  tympanal  membranes  (PTM,  see  Figure  1.A),  which 
present  a  relatively  large  surface  area  (compare  Figure  2.A  and 

2.B). 

2.  Inflation  of  the  ventral  probasisternum  (Pb)  to  provide  structural 
support  for  the  tympanal  membranes. 

3.  Bilateral  extensions  of  an  unpaired  sclerite  associated  with  the 
tympanal  membranes,  the  presternum  (Pr),  to  which  attach  the 
sensory  organs  at  the  tympanal  pits  (TP).  The  biomechanical 
function  of  the  presternum  is  a  key  to  the  process  of  directional 
hearing. 

4.  Enlargement  of  the  prosternal  air  sac  forming  the  acoustic  chamber 
backing  the  tympanal  membranes. 

5.  Location  of  the  two  scolopidial  sensory  organs  in  the  unpartitioned 
prosternal  air  sac.  A  scolopidium  is  the  elementary  multicellular 
arrangement,  including  a  mechanoreceptive  neuron  that  is 
commonly  used  in  insects  for  vibration  reception. 

6.  Cuticular  apodemes,  establishing  a  stiff  mechanical  link  between 
these  sensory  organs  and  the  presternum. 

7.  Reduction  in  size  of  the  prosternal  cervical  sclerites  (CSc). 

8.  Structural  reorganization  of  the  internal  endoskeleton  of  the 
prosternal  region  (Richard  R.  Fay  and  Arthur  N  Popper,  1998). 


6 


Figure  2.  Scanning  electron  micrographs  of  two  closely  related  species  of 
tachinid  flies  illustrating,  in  frontal  view,  the  anatomical  differences  of  the 
prosternal  region.  A.  The  tympanate  fly  Ormia  ochracea.  B.  The 
atympanate  fly  Myioppharus  doryphorae.  PeS,  proepisternal  setae:  PTM, 
prosternal  tympanal  membranes;  PM,  prosternal  membranes;  Pb, 
probasisternum;  Pr,  presternum;  Co,  prothoracic  coax;  N,  neck;  CSc, 
cervical  sclerite.  Notably,  in  the  tympanate  species  (A)  the  Pb  and  the  Pr 
are  conspicuously  larger,  and  the  PTMs  present  a  larger  surface  area  and 
are  thrown  with  radial  corrugations.  The  small  cuticular  depressions  at  the 
distal  ends  of  the  Pr  are  the  points  of  insertion  of  the  sensory  organs  with 
the  tympanal  system.  Scale  200  /urn  (Edgecomb  et  al.  1995;  Robert  et  al. 

1996a). 


A.  THE  MECHANICAL  ANALOGUE  OF  THE  INTERTYMPANAL 

COUPLING 

Analysis  of  the  deflection  shapes  of  the  tympanal  membranes,  and  more 
particularly  of  the  intertympanal  bridge,  indicates  that  the  two  ears,  join  at  the 
pivot  point.  The  incident  pressure  of  sound  produces  deflection  in  two  sides  of 
the  wing  shaped  tympana  two  deflections  of  wing,  which  are  almost  in  phase  and 
with  the  same  amplitude.  However,  as  shown  in  Figure  3,  the  intertympanal 
bridge  undergoes  an  asymmetrical  displacement  about  is  center,  much  like  a 
flexible  seesaw  rocking  back  and  forth  about  its  pivot  point. 

The  intertympanal  bridge  that  links  the  tympanal  membranes  is  a  key 
feature  of  the  auditory  mechanics  of  the  ormiine  ears.  The  ability  of  the 
intertympanal  bridge  to  rock  back  and  forth  in  a  flexible  manner  is  what  can  drive 
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both  membranes  with  such  phase  and  amplitude  differences.  The  ability  of  the 
incident  sound  pressures  to  drive  this  tympanal  system  depends  on  the  relative 
phase  of  the  pressures  acting  on  the  tympana.  In  the  same  way,  a  seesaw  is  put 
out  of  static  equilibrium  by  two  weights  of  equal  mass  applied  at  different  times 
on  each  of  its  ends;  these  forces  (pressures)  result  in  the  rocking  motion.  As  one 
arm  of  the  bridge  is  deflected  downward,  the  other  arm  will  move  upward  due  to 
the  stiffness  of  the  bridge  (Figure  3). 


Figure  3.  The  ears  of  O.  ochracea  and  a  mechanical  model  used  to  describe 
the  directional  sensitivity.  The  two  tympana  are  the  corrugated 
membranes  that  are  mechanically  connected  through  the  intertympanal 
bridge,  shown  here  with  the  numbers  1, 2,  and  3.  The  central  point  (3) 

acts  as  a  hinge. 


A  more  complete  analysis  of  this  system  shows  that  the  behavior  of  this 
mechanical  system  with  two  degrees  of  freedom  can  be  explained  by  the 
interaction  of  two  basic  modes  of  vibration  (one  rotational,  one  translational). 
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The  relative  contributions  of  these  modes  at  different  driving  frequencies  depend 
on  the  difference  or  sum  of  the  forces  acting  on  the  system  and  on  its  resonant 
properties  (see  Miles,  Robert,  and  Hoy  1995;  Robert,  Miles,  and  Hoy  1996b). 

This  morphological  and  biomechanical  evidence  led  to  the  development  of 
a  simple  mechanical  analogue,  the  mathematical  formulation  of  which  is 
presented  in  Miles,  Robert,  and  Hoy  (1995).  In  this  model,  the  flexible 
intertympanal  bridge  is  represented  by  two  rigid  bars  connected  medially  by  a 
torsional  spring,  and  dash  pot.  Both  extremities  of  the  bridge  are  connected  to  a 
spring  and  a  dash  pot  that  represent  the  stiffness  and  damping  characteristics  of 
the  auditory  sensory  organs  attached  to  the  tympanal  pits. 

The  deflection  shapes  computed  from  laser  vibrometer  and  COMSOL 
illustrate  the  unusual  rocking  and  transitional  motions  of  the  intertympanal  bridge. 
The  resulting  ipsilateral  and  contralateral  motion  amplitudes  can  be  visualized  in 
the  sketch  (Figure  4).  Unfortunately,  this  representation  does  not  illustrate  the 
phase  delay  introduced  by  this  floppy  connection. 


Floxiblo 

hingo 


N 


Figure  4.  Schematized  response  of  the  intertympanal  bridge. 


9 


However,  the  present  evidence  does  not  formally  exclude  that  an 
alternative  mechanism  could  account  for,  or  contribute  to,  the  observed  tympanal 
dynamics.  Indeed,  it  is  still  possible  that  the  air  space  behind  both  tympanal 
membranes  could  have  a  sufficient  stiffness  to  act  as  a  resonator  and  thus 
provide  the  basis  for  a  directional  response,  as  is  the  case  for  the  acoustically 
coupled  ears  of  pressure  difference  receivers  as  in  frogs  and  lizards.  In  a 
pressure  difference  system,  the  volume  of  the  air  chamber  backing  the  tympana 
is  crucial  for  the  directionality  to  be  determined  by  the  mechanical  response  of 
the  mechanical  response  of  the  tympana  (see  B.  Hedwig,  2005). 

B.  MULTI  DEGREE  OF  FREEDOM  VIBRATION 

In  this  section,  a  review  of  deterministic  vibration  is  provided.  These 
results  are  used  in  the  next  section,  where  the  forcing  is  taken  to  be  a  random 
process.  The  equations  of  motion  can  be  written  as: 

WW0}  +  W{x(0}  +  [^]M0}  =  W)}  (i.i) 

where  the  matrices  [m],  [c]  and  [k\  are  of  dimension  NxN ,  and  the  response 
{ x(t)}  and  force  {F(t)}  vectors  are  dimension  N  x  1 . 

This  concept  will  be  introduced  primarily  by  working  through  the  solution  of 
a  two  degree  of  freedom  system.  All  these  ideas  transfer  to  large  systems,  but 
with  the  two  degree  of  freedom  models,  they  can  be  demonstrated  without 
complications  of  the  major  algebraic  and  numerical  demands  made  by  the  larger 
systems. 
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Kt,  Ct 


Figure  5.  The  mechanical  model  includes  equivalent  stiffness,  Kt  and  Ks 

equivalent  viscous  dashpots,  C,  and  Cs . 
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Figure  6.  Free  body  diagram  of  a  two  Degree  of  Freedom  system  (see  Taylor 

&  Francis  Group  LLC,  2005). 

From  the  diagram  on  Figure  5  and  6,  motion  can  be  derived  using  the 
coupled  equation  of  motion  using  either  Newton’s  Second  Law  of  motion  applied 
to  a  free  body  diagram  for  each  mass  or  by  Lagrange’s  equation.  In  either  case, 
equations  reduce  to  governing  equations  to  be: 


mxl(t)  +  (Cj  +  c3)x1(0  +  c3X2(0+(^i  +  k1)xl(t)+  k3x2(t)  =  Fx(t)  (1 .2) 
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(1.3) 


mx2(t)  +  (c2  +  c3)x2(t)  +  c3xi(t)  +  k3x1(t)  +  (k2  +  k3)x2(t )  =  F2(t ) 


In  matrix  form  turn  to  be: 


ml 

0 


0 

m2 


x2(t)  J 


C,  +  C3 

C3 


^3 

r*2  ^3 


x2(t)  j 


k\  +  k3 
k3 


ki  ]k(Ol 
k2+  k3_  [x2(0j 


(1.4) 


For  this  case,  the  three  parameters  of  the  mass,  damping,  and  stiffness  are  given 
by 


0  ' 

\rl- 

c,  +  c3 

C3 

0 

111  t 

,  \c  J 

c3 

C2  +  C3  _ 

+  k^ 
k3 


k3 

h~  k ^ 


(1.5) 


C.  FREQUENCY  RESPONSE  FUNCTION 

Begin  by  taking  the  Fourier  transform  of  the  equations  of  motion,  obtaining 

{-a>2  [/«]+/'&>  [c]  +  [£j){X(<z>)}  =  (1.6) 

where  X(oo)  and  F(co)  are  Fourier  transforms  of  x(t)  and  F(t) .  The  matrix 
{-or  [m]+/ry[c]+  [£])  is  detonated  as  \Z{co  )].  Then 

{X{Q>)}=[Z{Q>)X  [F{(D)}  (1.7) 

The  matrix  [Z(®)] 1  is  identical  to  the  frequency  response  matrix  denoted  as 

[H{co)\ 
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[X(cd)}=[H{co)\{F(co)} 


(1.8) 


The  transfer  functions  are  given  by  H np(co)  =  seim/1  and  Hf2p(co)  =  se  icox/1 , 

where  z  is  the  time  for  the  incident  sound  to  travel  between  the  points  where  the 
forces  act  as  given  above  and  s  is  the  surface  area  of  each  tympanal  membrane 
as  shown  below. 


Figure  7.  Time  delay  between  the  ipsilateral  and  contralateral  mechanical 
response,  calculated  from  the  difference  phase  spectrum  from  sound 

source  at  incident  angle. 


The  distance  of  d  is  the  distance  sound  travel  with  respects  to  the  angle  6 
and  c  (speed  of  sound).  The  time  delay  of  r  can  be  determined  as: 


d  Lsm6 
c  c 


(1.9) 


The  transfer  functions  between  the  responses  and  the  pressure  at  the 
pivot  may  be  found  using  Cramer’s  rules,  a  method  to  solve  systems  equations 
using  determinants.  Starting  with  equation  below: 


2 

m ! 

0  ' 

ct+c3 

C3 

X1  , 

kx  +  k3 

k3 

xl 

Fi 

CO 

0 

m2_ 

+  ICO 

x2 

C3 

C2  +  C3  _ 

+ 

x2 

1 

UJ 

k2  +  k3  _ 

x2 

ii 
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(1.11) 


-co2ml  -  ico(cl  +  c3)  +  kx  +  k3 


- icoc 3  +  k3 


- icoc 3  +  k3 


-co2m2  -  ico(c2  +  c3)  +  k2  +  k3 


F, 

F, 


By  Cramer’s  rules  method,  the  transfer  function  between  H 
FfXiP(co )  the  responses  and  the  pressure  at  the  pivot  are  shown  below: 


det 


HXlP(co)  = 


F 


-icoc3  +  L 


F2  -co2m2  -  ico(c2  +  c3)  +  k2  +  k3 


det 


-co2ml  -  ico(cl  +  c3)  +  kx  +  k3 


-icoc3  +  k3 


- icoc3  +  k3 


-co1  m2  -  ico{c2  +  c3)  +  k2  +  k3 


s{k3  +  icoc3)x  (e  -e  )  +  s(k  +  icoc- mco~)e 

Hx  ( CO )  = - - - — - ; — 

(k  +  icoc  +  k3  +  icoc3- mco  )  -{k3  +  icoc3Y 


2  \  icozll 


det 


HXlP{co)  = 


-co2m1  -  ico(C\  +  c3)  +  kx  +  k3  Fl 


-icoc3  +  k3 


F 


det 


-icoc3  +  k3 

- icoc „  +  k ,  -co2 m2  -  ico{c2  +  c3)  +  k2  +  k3 

2  \  -icor/2 


-co  ml  -  ico{cx  +  c3)  +  kx  +  k3 

'3  +  ^3 

-icozt  2  -icQT/2 


HXlP{oo)  = 


s{k3  +  icoc3 )  x  ( e  U0T  “  -  e  U0T  )  +  s{k  +  icoc  -  mco  )e 
(k  +  icoc  +  k3+  icoc3  -  mco2)1  -  ( k3  +  icoc3 )2 


where, 


Rocking  frequency,  cor 


Bending  frequency,  co, 


(k  +  2  k3) 


m 


{co)  and 

(1.12) 

(1.13) 

(1.14) 

(1.15) 

(1.16) 

(1.17) 
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c 


(1.18) 


£  = 


2  corm  ’ 


(c  +  2c3) 
2o)fm 


To  simplify  the  notation,  let  Hfp(co)  and  Hfp(a> )  denote  the  transfer 

functions  as  a  function  of  the  frequency,  co,  in  radians/s  of  the  ipsilateral  and 
contralateral  force  fx(t),  and  f2(t) ,  relative  to  the  acoustic  pressure  at  the 
microphone  location,  p{t).  The  solutions  of  modes  for  the  transfer  function  in 
Equation. 1 .12  and  1.14  are  given  by: 


HXlP(co)  = 


(IIfiP(co)-IIfJco))/m 
co2.  -  co1  +  4corgrico 


+ 


(fIfiP(co)+  fI/iP(  CO))1  in 
cof  -  co2  +  4cot^ico 


(1.19) 


HX2P{(o)  = 


(flfP(co)+  nf2p(  co))1  in 
co 2  -  co2  +  4 co^ico 


co2  -  co2  +  4 calico 


(1.20) 


where  Hxp(co)  and  Hx  p(co)  are  the  transfer  functions  between  the  responses  of 
the  ends  of  the  tympanal  bridge,  xx(t )  and  x2(t) ,  relative  to  the  pressure  at  the 
pivot  point,  p(t)  (see  Miles,  Robert,  and  Hoy  1995). 

D.  THE  EAR’S  SENSITIVITY  TO  THE  DIRECTION  OF  SOUND 

The  expression  of  the  two  oscillation  are  described  as  in  Equation  1 .20  to 
1.22  which  is  explained  at  Figure  5  and  8  for  two  Eigen  modes  (G.  Karunasiri, 
2008).  The  displacements  of  the  two  sides  can  be  expressed  as: 


.  ( <z>dsin<9^ 

.  ( <Z)dsin^ 

Ps 

sm 

l  2VS  J 

A 

sm 

l  2VS  J 
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yl(co2-co2J  +(rTco)2 
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Af  =  A;  +  A;  +  2 AbAr  sin(^b  -  £),  A,  =  Ab  +  A;2  -  2AbAr  sin(^b  -  fc) 


(1.22) 


Xj  =  A,  sin(<yt  +  <j>x)  ,  x2  =  A2  sin(<yt  +  $,)  (1 .23) 

where  A,  and  A2  are  the  amplitudes  of  ipsilateral  and  contralateral  and  ^  and  ^2 
are  the  corresponding  phases  (Figure  8).  The  sound  direction  can  be 
determined  by  looking  at  the  amplitude  and  the  phase  difference  between 
ipsilateral  and  contralateral. 


Figure  8.  Phase  difference  and  ratio  of  amplitude  are  independent  of  sound 

pressure. 
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III.  SENSOR  DESIGN 


The  4th  generation  sensor  chip  sensor  chip  contains  12  different  sensors 
as  shown  in  Figure  9.  The  chip  was  fabricated  using  SOIMUMPs  process 
available  through  the  MEMSCAP  foundry  service  (SOIMUMPs  Design 
Handbook,  Rev.  4.0). 


Figure  9.  Layout  of  the  4th  generation  chip.  The  Devices  4  and  6  are 

characterized  in  this  work. 
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In  the  following,  details  of  the  Devices  4  and  6  in  the  layout  shown  in 
Figure  10  are  given  since  they  are  extensively  characterized  in  the  work. 
Dimensionally  they  are  identical  except  the  Device  4  the  wings  are  perforated 
with  square  holes  having  of  2  x  2  |um2  area.  The  Table  1  gives  the  dimensions  of 
the  structure. 


E 

& 

A 

C  C02 

B 

I 

D 

R1 

R2 

Figure  10.  Sensor  is  designed  to  achieve  the  desirable  operating  frequency. 


A  and  B 

C 

D  and  E 

Dimension 

1  mm  x  0.75  mm 

0.080  mm  x  0.5  mm 

0.045  mm  x  0.075  mm 

Table  1 .  The  sensor  dimension  for  both  devices. 


A.  FABRICATED  SENSOR  CHIP  (DESIGNED  USING  MEMSPRO) 

Figure  1 1  shows  the  expanded  view  of  the  two  devices.  The  substrate 
under  the  devices  was  trenched  during  the  processing  to  reduce  squeezed  film 
damping  as  depicted  in  Figure  12. 
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(a) 


(b) 


Figure  1 1 .  The  characteristics  of  both  wings  of  each  device.  A.  Device  4 
(perforated  wings).  B.  Device  6  (solid  wings). 


Figure  12.  The  air  gaps  between  the  silicon  substrate  and  the  wings. 

One  of  the  main  goals  of  the  research  is  the  probe  the  effects  of  holes  on 
the  device  4  on  the  dynamics  of  the  sensor  response.  In  particular,  the  coupling 
of  sound  and  enhanced  damping  due  to  possible  airflow  through  the  holes. 

B.  EXPERIMENT  SETUP 

Figure  13  shows  a  picture  of  the  laboratory  setup  used  to  measure  the 
vibrational  amplitudes  of  the  sensors  under  sound  excitation.  The  experiment 
setup  consists  of  the  reference  microphone,  speaker  on  a  rotating  boom,  laser 
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vibrometer,  and  devices.  The  reference  microphone  was  placed  close  to  the 
sensor  chip  to  measure  the  sound  pressure.  It  is  a  Bruel  &  Kjaer  Pressure  field 

—  microphone  type  4138.  It  has  a  relatively  flat  response  curve  from  20  Hz  to  20 

8 

kHz  with  0.939  mV/Pa  sensitivity. 

The  sound  source  was  Selenium  loudspeaker  type  DH200E  attached  to 
the  internal  signal  generator  in  the  VibSoft  software.  The  software  allows  the 
generation  of  pure  tones  as  well  as  chirp  sound  with  different  amplitude  and 
frequency.  Incident  sound  wave  angle  was  adjusted  manually  by  adjusting  the 
speaker  position.  The  laser  vibrometer  was  a  Polytec  single  point  vibrometer 
model  OFV  302,  with  a  model  OFV  2600  controller.  Its  purpose  was  to  measure 
the  displacement  of  the  wings,  to  a  resolution  in  the  sub  nanometer  range. 
(Antonios  Dritsas,  thesis  paper,  June  2008). 


Figure  13.  The  existing  equipment  used  for  measuring  vibrational  amplitudes. 
A.  Rotating  sound  source.  B.  Sensor  devices  and  microphone.  C.  Speaker 
for  sound  source  with  directional  cone.  D.  Laser  vibrometer. 
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The  sound  source  attached  to  the  rotating  stage  in  Figure  14  is  used  to 
measure  the  response  at  a  discrete  set  of  angles. 


Figure  14.  Rotating  stage  for  the  sound  source. 


The  location  of  the  reference  microphone  relative  to  the  sensor  chip  is 
shown  in  Figure  15.  The  microphone  is  connected  to  the  reference  input  of  the 
laser  vibrometer. 


Figure  15.  Sensor  and  microphone. 


The  speaker  (see  Figure  16)  used  as  the  sound  source  has  a  relatively  flat 
response  in  the  1-2  kHz  frequency  rage  consists  of  a  cone  to  direct  the  sound 
towards  the  sensor. 


21 


Figure  16.  Speaker  for  sound  source  with  directional  cone. 

The  laser  vibrometer  shown  in  Figure  17  is  to  measure  the  vibration 
amplitudes  of  the  wings  in  response  to  incident  sound  at  different  angles.  The 
typical  displacements  of  the  wings  are  in  the  tens  of  nanometers. 


Figure  17.  Laser  vibrometer  for  measuring  displacement  amplitudes  under 

sound  excitation. 


C.  MEASURED  FREQUENCY  RESPONSE  (DEVICE  4  AND  6) 

The  deflection  amplitudes  computed  from  the  laser  vibrometer 
measurements  illustrate  the  expected  rocking  and  bending  motions  of  the  wings 
of  the  sensors.  The  measured  rocking  and  bending  motion  amplitudes  as  a 
function  of  frequency  for  the  two  devices  at  45°  are  shown  in  the  Figure  18.  The 
measured  amplitude  at  the  bending  mode  is  about  10  pm  per  1  Pa  of  sound 
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pressure.  The  amplitudes  at  the  rocking  motion  were  relatively  small  since  it 
depends  on  the  arrival  time  difference  between  the  two  wings,  which  is  order  of  a 
few  microseconds. 


Figure  18.  The  measured  frequency  response  for  Device  4  (perforated  wings) 

and  6  (solid  wings). 

The  sound  intensity  is  set  by  applying  0.1  Volt  to  the  controller,  which  is 
connected  to  an  amplifier  for  driving  the  speaker.  The  amplifier  generates  an 
output  voltage  depending  on  the  gain  used  to  drive  the  sound  source.  The  signal 
from  the  reference  microphone  was  converted  to  pressure  using  its  conversion 
factor  of  0.939  V/Pa.  The  measured  rocking  frequency  of  the  Device  4  and 
Device  6  were  approximately  at  4385  Hz  and  3680  Hz,  respectively  while  the 
bending  mode  are  found  to  be  around  4850  Hz  and  4620  Hz  as  shown  in  Figure 
18.  The  two  devices  show  different  resonant  frequencies  as  well  as  peak 
amplitude.  The  Device  4  (perforated  wings)  shows  slightly  smaller  amplitude  of 
deflection  for  both  the  rocking  and  bending  modes  compared  to  that  of  the 
Device  6  (solid  wings)  as  depicted  in  Figure  20. 
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In  experiment,  the  response  of  the  sensors  was  also  investigated  by 
random  noise  analysis.  In  this  analysis,  the  stimulus  was  a  band  limited  random 
noise  burst  (500  to  8000  Hz  bandwidth)  delivered  at  45°  of  incidence  angle.  The 
mechanical  responses  of  the  ipsilateral  and  contralateral  are  measured  by 
positioning  the  beam  of  the  laser  vibrometer  on  the  locations  indicated  in  Figure 
3.  The  ipsilateral  side  found  to  vibrate  in  the  sound  field  with  slightly  higher 
amplitude  than  the  contralateral  side. 

D.  SIMULATED  FREQUENCY  RESPONSE  (DEVICE  4  AND  6) 

In  order  to  compare  the  measured  frequency  responses  with  the  designed 
values,  a  finite  element  analysis  was  carried  out  using  COMSOL.  Figure  19 
shows  the  simulated  rocking  and  bending  motions  of  the  sensor  with  solid  wings. 


Figure  19.  Deflection  shapes  of  (a)  rocking  and  (b)  bending  modes. 


The  simulated  frequency  responses  of  the  two  devices  show  different 
peak  amplitudes  and  frequencies.  Device  4  (perforated  wings)  shows  smaller 
deflection  than  Device  6  (solid  wings)  in  both  modes  as  shown  in  Figure  20.  The 
simulated  frequency  responses  are  in  good  agreement  with  the  measured  data 
shown  in  Figure  18. 
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Figure  20.  The  simulated  frequency  response  for  Device  4  (perforated  wings) 

and  6  (solid  wings). 
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IV.  THE  MECHANICAL  TO  ELECTRICAL  TRANSFORMATION 

(PSPICE) 


In  order  to  determine  the  speed  of  operation  of  the  sensors,  the  transient 
response  of  them  were  analyzed  using  SPICE  circuit  simulator  based  on 
electrical  and  mechanical  equivalents.  The  following  briefly  discusses  the 
analogies  between  electrical  and  mechanical  components.  Table  2  lists  the 
correspondence  between  mechanical  and  electrical  quantities. 


MECHANICAL  QUANTITY 

ELECTRICAL  QUANTITY 

Force,  F 

Electromotive  force,  V 

Velocity,  u 

Current,  / 

Displacement,  x  =  J udt 

Charge,  q  =  Jm 

Mass,  M 

Inductance,  L 

Mechanical  Resistance,  Rm 

Resistance,  R 

Stiffness,  km 

Capacitance,  1/C 

Table  2.  Electrical-Mechanical  equivalencies. 


Based  on  the  equivalent  quantities  shown  in  Table  2,  the  electrical 
equivalent  circuit  of  the  mechanical  model  of  the  fly's  hearing  system  in  Figure  5 
can  be  constructed  as  shown  in  Figure  21 . 


LI  R1  Cl  C2  R2  L2 


Figure  21 .  Electrical  equivalent  circuit  of  the  mechanical  model  of  the  fly's 

hearing  system  in  Figure  5. 
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The  differential  equations  of  motion  for  the  mechanical  system  and  the 
electrical  circuit  have  the  same  form  as  shown  by  Equations  1 .24  -  1 .27. 


Mij+  (Cj  +  C3)Xj+  C3x2  +  (£,  +  k2)xl  +  k2x2  =  F 


(1.24) 


d2 ql 


dqx  ,  D  dq2  T  ,  f  1 , 0  „  ,  1 


R~t^~  +  (Ri  +  Ri)~  +  +  +  +  ~q2  ~  v\ 


dt 


dt  dt 


C, 


(1.25) 


M  x2+  C3  X[+  (C2  +  C3)x2+  A:3X[  +  (Z2  +  Z3)x,  =  F2 


(1.26) 


d2  q0  „  dq,  ,  „  „  .  dq ,  1  ft  1  ^ 

Z - r — I-  i?3 - 1-  (i?„  +  i?3) — —  H - qx  +  I - 1 - I  q2  = 

dt 2  3  dt  2  3  dt  c3  1  Ic,  cj  2  2 


(1.27) 


The  analogy  is  very  useful  because  the  equivalent  electrical  circuit  can  be 
readily  solved  using  PSPICE.  The  parameter  of  the  electrical  circuit  was 
obtained  from  the  measured  frequency  responses  of  the  two  sensors  given  in 
Figure  21.  The  mass  of  each  wing  was  estimated  based  on  the  dimension  and 
density  of  Si  and  the  mechanical  resistance  was  estimated  using  the  peak  widths 
and  Tables  3  and  4  list  the  parameters  for  the  sensors  with  perforated  and  solid 
wings,  respectively. 
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A.  IMPEDANCE  ANOLOGY  FOR  DEVICE  4  (PERFORATED  WINGS) 


Mechanical  Quantity 

Electrical  Quantity 

Remarks 

Mechanical 

Resistance 

Rrrn  and  Rrri2 

( 3  x  1 0  6  kg/s  ) 

Resistor 

Ri and  R2 
( 3  x  10~6  Q) 

Rmi  and  Rrri2  = 

Ri and  R2 

Rrri3 

( 5.2 x  1CT6  kg/s) 

r3 

(5.2x10  6  Q) 

Rm3  =  R3 

Mass 

Mi  and  M2 

(4.66 x  1CT9  kg/s) 

Inductor 

Li  and  L2 
( 4.66  x  10~9  Q) 

Mi  and  M2  = 

Li  and  L2 

Stiffness 

ki  and  k2 

(3.71  N/m) 

Capacitance 

Ci and  C2 
(0.27  F) 

ki  and  k2  = 

—  and  — 

C  C 

k3 

( 0.685  N/m ) 

Cs 

(1.46  F) 

k,  =  — 

c 

v-'3 

Table  3.  Mechanical  and  electrical  quantities  for  Device  4. 


B.  IMPEDANCE  ANOLOGY  FOR  DEVICE  6  (SOLID  WINGS) 


Mechanical  Quantity 

Electrical  Quantity 

Remarks 

Mechanical 

Resistance 

Rmi  and  Rrrn 

(2.93x  10  6  kg/s) 

Resistor 

Ri and  R2 

(2.93x  10  6  Q ) 

Rmi  and  Rrri2  = 

Ri and  R2 

Rrri3 

( 4.2 x  10  6  kg/s) 

r3 

(4.2x10  6  Q) 

Rm3  =  R3 

Mass 

Mi  and  M2 

(4.66 x  10~9  kg/s) 

Inductor 

Li  and  L2 
(4.66  x  10  9  Q) 

Mi  and  M2  = 

Li  and  L2 

Stiffness 

ki  and  k2 

(2.632  N/m) 

Capacitance 

Ci and  C2 
(0.38  F) 

ki  and  k2  = 

—  and  — 

c  c 

k3 

( 0.526  N/m ) 

C3 

(1.9  F) 

k  =  — 

c3 

Table  4.  Mechanical  and  electrical  quantities  for  Device  6. 
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This  graph  in  Figure  22  shows  the  calculated  amplitudes  using  PSPICE 
based  on  the  circuit  described  in  Figure  21.  The  amplitudes  were  estimated 
based  on  the  charge  accumulated  on  the  capacitors.  For  the  Device  4,  the 
rocking  amplitude  is  about  0.76  pm  (CXV  =  021V  x  2.8//V)  and  12.41  pm  for  the 
bending  frequency,  (C3F=1.46F  x  8.5//V).  The  corresponding  values  for  the 
Device  6  were  found  to  be  about  1 .06  pm  for  rocking  frequency  and  10.5  pm  for 
bending  frequency. 


Figure  22.  Simulated  frequency  response  at  45°  using  PSPICE. 

C.  PEAK  HEIGHT,  PEAK  WIDTH  AND  QUALITY  FACTOR 

The  COMSOL  simulated  and  measured  amplitudes,  peak  widths  and 
quality  factors  (Q)  for  the  Devices  4  (perforated  wings)  and  6  (solid  wings)  are 
shown  in  Table  5.  The  width  of  the  peak  is  defined  as  the  distance  between  the 
two  half-power  points  at  +  3dB frequencies  as  in  Figure  23.  The  distance 
between  half  power  points  is  referred  to  as  the  bandwidth: 
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Peak  width  =  /+3dB-/_3dB 


(1.28) 


The  quality  factor  is  a  term  used  to  describe  the  sharpness  of  the  peak  as  it 
defines  by: 


Quality  Factor, 


0  = 


fo 


/+3dB  -/-  3dB 


(1.29) 


Figure  23.  The  peak  width  of  rocking  frequency. 


COMSOL 

LASER  VIBROMETER 

Device  4 

Device  6 

Device  4 

Device  6 

Peak  height 

0.707  pm 

0.742  pm 

0.636  pm 

0.606  pm 

Peak  width 

34  Hz 

23  Hz 

35  Hz 

17  Hz 

Quality  Factor 

131 

168 

123 

214 

Table  5.  The  comparison  of  two  devices  in  term  of  peak  height  width,  and 

quality  factor. 
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D. 


TRANSIENT  TIME  OF  RESPONSE 


The  transient  responses  of  the  two  sensors  were  simulated  using  PSPICE 
the  sound  incident  at  45°.  Figures  24  (a)  and  (b)  show  the  simulated  transient 
data  for  the  Device  6  and  4,  respectively. 


Figure  24.  Transient  time  of  response  measured  by  PSPICE.  A.  Device  6 
(solid  wings).  B.  Device  4  (perforated  wings). 

The  data  in  Figure  24  show  that  the  time  required  for  attaining  the  steady 
state  response  for  Device  4  (perforated  wings)  is  about  2.2  ms  and  for  Device  6 
(solid  plate)  is  about  3.1  ms.  The  higher  settling  time  for  the  Device  4  is  due  to 
the  smaller  damping  ratio  y  =  b/m.  The  transient  response  amplitude  of  a 
damped  harmonic  oscillator  is  given  by  (G.  Karunasiri,  2009) 

A(t)  =  A0eyt  (1.31) 
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The  lighter  mass  associated  with  the  structure  with  perforated  wings 
produce  a  larger  damping  ratio  making  a  faster  rise  of  the  transient  amplitude.  In 
addition,  the  perforation  can  increase  the  damping  coefficient  (b)  making  the 
damping  ration  higher. 

E.  ROCKING  AND  BENDING  AMPLITUDE  RATIO  WITH  INCIDENCE 
ANGLE  OF  SOUND 

In  order  to  determine  the  directional  sensing  capability  of  this  system, 
amplitudes  at  rocking  and  bending  frequencies  were  measured  for  as  a  function 
of  direction  of  sound.  Figure  25  shows  the  simulated  dependence  of  ratio  of  the 
rocking  to  bending  amplitudes  using  COMSOL  with  the  incident  angle  for  the  two 
sensors.  As  expected,  the  ratio  increases  with  angle  due  to  the  increase  of  the 
arrival  time  delay,  which  increases  the  rocking  amplitude.  The  rationing  was 
used  to  remove  the  unknown  sound  pressure  at  the  sensor  since  the  increase  in 
time  delay  does  not  affect  the  bending  motion. 


Figure  25.  Directional  sensitivity  of  the  mechanical  response  from  COMSOL. 
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The  Table  6  shows  the  measured  response  for  Device  4  (perforated 
wings)  and  Device  6  (solid  wings),  using  COMSOL  for  a  number  of  different 
incident  angles  of  sound.  All  the  measurements  were  taken  with  intensity  of 
sound  wave  of  0.939  mV/Pa.  The  increase  of  the  ration  with  angle  agrees  well 
with  the  simulated  data  shown  in  Figure  25. 


Phase 

(°) 

DEVICE  4  (perforated  wings) 

DEVICE  6  (solid  wings) 

Amplitude  (pm) 

Ratio 

Amplitude  (pm) 

Ratio 

Rocking 

Bending 

Rocking/Bending 

Rocking 

Bending 

Rocking/Bending 

10 

0.43 

3.7 

0.116 

0.38 

8.15 

0.0466 

20 

0.52 

3.7 

0.141 

0.44 

8.15 

0.0539 

30 

0.60 

3.7 

0.162 

0.49 

8.15 

0.0601 

40 

0.68 

3.7 

0.183 

0.53 

8.15 

0.0651 

50 

0.75 

3.7 

0.202 

0.57 

8.15 

0.0699 

60 

0.79 

3.7 

0.214 

0.60 

8.15 

0.0736 

70 

0.83 

3.7 

0.224 

0.63 

8.15 

0.0773 

80 

0.87 

3.7 

0.235 

0.64 

8.15 

0.0785 

Table  6.  Difference  in  the  amplitude  of  the  mechanical  response  between 

the  tympanal  pits  for  different  angles  of  sound  incidence. 
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V.  CONCLUSION 


Based  on  the  structure  of  the  auditory  system  of  the  fly  Ormia  ochracea,  a 
directional  microphone  was  designed  and  characterized.  It  was  designed  to 
fabricate  using  SOIMUMPs  the  performance  was  analyzed  using  COMSOL  and 
PSIPCE  software.  The  main  goal  of  the  research  was  to  characterize  two 
identical  sensors  with  solid  and  perforated  wings.  It  was  found  that  the  two 
sensor  had  nearly  identical  response  amplitude  to  sound  indicating  good  sound 
coupling.  The  sensor  with  perforated  wings  showed  less  curling  due  to  residual 
stress  compared  to  the  sensor  with  solid  wings.  Both  sensors  showed  good 
sensitivity  to  the  direction  of  sound.  The  measured  and  simulated  frequency 
responses  of  the  sensors  showed  good  agreement.  The  use  of  perforated  wing 
allows  the  increase  of  damping  ratio  making  the  sensor  to  respond  faster  and 
giving  a  relatively  broader  frequency  response. 
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VI.  RECOMMENDATIONS  FOR  FUTURE  WORK 


For  the  future  work,  these  experiments  need  to  be  done  in  an  anechoic 
chamber  to  reduce  the  interference  from  unwanted  reflection  from  the  walls  in  the 
lab. 
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